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Airfoil chord length, m 
Drag coefficient, D/q^c 
Lift coefficient, L/q^c 

Moment coefficient about the quarter chord, M/q c2 
Pressure coefficient, (P - P 00 )/q 00 
Liquid water content, g/m3 
Local static pressure, N/m^ 

Free stream static pressure, N/m 2 

p 

Free stream dynamic pressure, N/m 
Temperature, °C 
Time, minutes 

Free stream velocity, m/sec 
Angle of attack, degrees 
Particle diameter, microns 
Air density, Kg/m^ 

Ice density, g/cm^ 
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INTRODUCTION 


Until recently very little experimental data has been 
available on the performance degradation of airfoil sections 
resulting from rime or glaze ice accretions. The test program 
described in this report provides some of this needed infonnation 
It's primary objectives were: 

1) To examine a method of simulating ice accretions with wood 
shapes which were instrumented with surface pressure taps 
to obtain detailed aerodynamic data. 

2) To study the complex flowfield in the area of the ice 
accretion through these pressure distributions. 

3) To evaluate the accuracy of the theoretical analysis methods 
presently utilized 1 . 

This paper is primarily a data report and presents informa- 
tion taken on the 63A**15 airfoil. Further analysis on other 
classes of airfoils is required to fully analyze the icing prob- 
lem and to refine conclusions drawn from the results of this test 

Background 

The first tunnel entry, conducted May **-15, 1981, was an 
ice accretion study. Rime and glaze ice shapes were grown on a 
NACA 632”A**15 airfoil with a 1.36 m chord. The resultant section 
drag coefficients were measured using a drag wake survey probe 
positioned approximately one chord length downstream of the model 
Tests were conducted at free stream conditions to simulate 
typical aircraft cruise (high velocity, low angle of attack) and 
climb (low velocity, high angle of attack) conditions. In order 
to generate glaze Ice a total temperature of -**°C was chosen 
while -26°C was used to yield rime st ipes. 
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Two methods were used to record the characteristics of the 
accreted shapes. The first involved fitting a template to the 
leading edge of the wing after a small section of the ice near 
the centerline had been removed. A tracing of the shape could 
then be made onto the template. The second method, typically 
used for longer icing times, involved removing a section of the 
accretion by spraying steam inside the model near the area to be 
removed. The ice section could then be lifted out and dipped 
into a container of beeswax which hardened around the ice. After 
the ice has melted, plaster was poured inside and the casts were 

p 

then available for careful tracings . 

Prom the shapes generated in the first entry, 2 rime and 
2 glaze shapes were chosen to represent typical climb and cruise 
conditions. These were then modelled for the second tunnel entry. 
Table I gives a summary of the test parameters which produced 
these shapes. 


TABLE I 

TEST PARAMETERS 


TYPE 

T 

a 

U 

6 

LWC 

t 

p ice 

Rime 

-26 

2.6 

51 

15 

1.5 

15 

0.421 

Rime 

-26 

6.6 

40 

15 

1.5 

15 

0.53*J 

Glaze 

- 4 

2.6 

51 

15 

1.5 

15 

- 

Glaze 

- 4 

6.6 

40 

20 

2.9 

15 

- 
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Equipment 

The testing was performed In the NASA Lewis 6’ x 9’ Icing 
Research Tunnel. Tracings of the shapes generated during the 
first tunnel entry were used to make the simulated ice shapes 
for the second test. The shapes were formed from mahogany and 
extended full span with filletting at the ends for attachment 
to the model. The Inside of each of the four shapes, two rime 
and two glaze, was hollowed out to provide clearance for the 1/16" 
ID tubing required for tapping the ice shape. In an effort to 
obtain as much data as possible in the critical region around the 
leading edge, the number of taps on the shape was concentrated in 
this area (see Figures 1A-1D). The airfoil was instrumented to 
obtain surface pressure measurements, using 1/8" OD strip-a-tube 
attached to the airfoil surface (Figure 2). In order to add sur- 
face roughness to the ice shapes, aluminum oxide grit with a 
K/C = .00058 was attached using a spray acrylic adhesive to the 
glaze shapes while for the rime cases a grit with K/C = .0012 
was used. 

Data acquisition and on-line reduction was accomplished 
using the Digital Data Acquisition and Reduction System (DDARS)3, 
designed at The Ohio State University’s Aeronautic il and Astro- 
nautical Research Laboratory. The DDARS package is illustrated 
in Figure 3 . The heart of the system is the central px’ocessing 
unit (CPU), which is a Digital Equipment Corporation (DEC) LSI-11 
microcomputer. The analog data signals from the various trans- 
ducers and the wake probe slldewire are routed to the computer 
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through the analog front-end, consisting of the various analog 
and signal conditioning equipment. Input and output is accomp- 
lished through a teletype terminal and data storage is achieved 
through a dual-drive single head double-density floppy disc 
system. 

In order to obtain pressures on the airfoil surface a 
Scanivalve transducer set-up was used (Figure 4). Connected 
through an analog box to the LSI- 11 microcomputer were wake 
toval, static, tunnel velocity, and tunnel static transducers 
taps. Drag data were taken using the wake momentum deficit 
technique. While the drag data obtained during the first tunnel 
entry were obtained using the IRT wake probe, for the aerodynamic 
testing an airfoil mounted probe was used. It was attached to 
the lodel on centerline and 0.413m back from the trailing edge, 
and measured both wake total and wake static pressure. 

Data Reduction 

On-line reduction is a key feature of the OSU DDARS. A 
quick-look C p distribution is made available to the engineer 
as well as integrated values of C£, Cm, and C(j. This permits 
immediate evaluation of the progress of the test and thus maxi- 
mum use can be made of tunnel time. 

Final reduction was performed on the OSU Harrls/6 Computer 
System. The raw data obtained daring testing in the IRT were 
copied onto the Harris System from the LSI-11. Pressure coef- 
ficient distributions, as well as lift and moment coefficients 
were obtainted In the same manner as the on-line reduction, how- 
ever a provision in the program permitted omitting from t-re 
final results any blocked taps on the model. 
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The drag coefficient reduction program displayed the wake 
on a graphics t-.erminal. This allows the operator to choose the 
integration limits required by the Jones Equation^. If it was 
clear from the displayed wake trace that the extent of the 
probe’s traverse was not sufficient to capture the entire wake, 
that reduction was bypassed. This limitation was observed at 
the higher angles of attack and is represented in the run summary 
tables by a dashed line in the drag coefficient column. 
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During the course of the wind tunnel test program a few 
observations were made which merit discussion and are pertinent 
to further operation of the Icing Research Tunnel. These com- 
ments will be made first. Following, a brief summary of the 
data collected is presented. The tabulated Rum Summary and 
plots of Cp’s and integrated results are given in the appendices. 

IRT Operation 

The initial aerodynamic tests were conducted with the model 
installed in the wind tunnel with a 0.06m (2 inch) gap between 
the end of the wing and the tunnel ceiling. This gap, at high 
lift coefficients, can influence the quality of the desired two- 
dimensional data; therefore, any clearances between two-dimen- 
sional airfoil model ends and the tunnel walls must be minimized. 

The early tests (runs 2 to 13) with the gap present can be 
compared with test conducted on the model. Figure 2, with sealed 
ends in Figure 5, 6 and 7* These figures present the lift coef- 
ficient vs. angle of attack, drag polar, and moment coefficient 
vs. lift, respectively at Reynolds number near 5 x 1U^. The 
solid line shows on the figures represents earlier NACA data of 
the 632”^15 airfoil at Reynolds number of 6 x 10^ with NACA 
standard roughness applied to the leading edge. 

Data scatter in the vs. a plot of Figure 5 is attributed 
to difficult; angle of attack measurements, but at high 

angles the lift coeff. cients of the gap tests appear below those 
of the sealed gap. Unfortunately, angles above 11° were not ex- 
amined for the gap case. Figure 6, the drag polar, removes the 
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angle of attack discrepancies and clearly illustrates the devia- 
tions between the gap and sealed gap tests. The open symbols, 

representing the sealed gap, track the early two dimensional NACA 
5 

data , while the gap data at high lift appear to have low drag 
coefficients. The anomalous behavior is caused by the relative 
locations of the pressure belt used to determine lift coefficient 
and the wake probe that measured the drag coefficient. The prob*> 
was near the tunnel centerline, while the belt was inboard 0.2 m 
further from the gap. The effect of the gap is to alter the 
spanwise angle of attack, producing a lift coefficient at the 
belt location greater than the lift coefficient at the probe 
station, and in turn generating the drag polar of Figure 6. A 
truly two dimensional test would, of course, have a constant 
spanwise angle of attack. The moment coefficient curves of 
Figure 7 illustrate a similar disagreement at high lift coef- 
ficients . 

A second observation made during the test program was the 
need for more accurate angle of attack setting and monitoring. 

The scatter in the vs. a plot is largely due to the inability 
to set the angle ef* the turntable to better than +0.1 degree. 

A* a result of these comparisons, the need to ensure that 
the model spans the entire tunnel has been demonstrated. Only 
the data obtained with the gap sealed will be used for later 
comparisons. Also, the problem with the tunnel turntable in 
setting and measuring the angle of attack has been discussed. 

The purpose of the test program was in part to evaluate and 
demonstrate *he use of the IRT as an aerodynamic tunnel . The 
problem with the turntable has been identified and can hopefully 
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be Improved. Other possible problems such as tunnel flow quality 
are still being studied. It is important then to note when using 
these data, that the IRT is still being improved, and is under 
evaluation as an aerodynamic facility. 

Simulated Ice Result s 

Comparisons wer a made between the drag measured on the 
actual iced airfoil and the simulated shapes in order to deter- 
mine the accuracy of the model. In addition, experimental 
values of C i and Cm as well as pressure distributions were com- 
pared with theoretical results^. 

Data were taken on the four smooth and rough simulated 
configurations, as well as the clean airfoil for lift coefficients 
of zero up to The glaze ice C p distributions show very 

clearly the adverse pressure gradient in the region around the 
leading edge. Prominent pressure spikes can be seen on the dis- 
tribution where the flow is forced around the horns of the ice 
shapes. These severe adverse gradients will promote separation 
and lower the C£ max , as well as increasing the drag of the wing. 
This problem was also seen in the rime cases, however, it was 
not as severe. 

Interestingly enough, even though the simulated shapes do 
show a reduction in C£ max , it was not as large as expected. In 
general, the rime cases failed to show a reduction and actually 
increased C£ max by 0.2 in the climb case, essentially forming a 
leading edge flap. The glaze shapes reduced C£ max by approxi- 
mately 0.2. 

The drag results for the Rime 3 case show the importance 
of the addition of the surface roughness. The comparison between 


PAGE 13 
QUALITY 


ft 


ORIGINAL PAGE 13 
OF POOR QUALITY 

the actual ice and simulated Ice was very good when the roughness 
was added. However, the Rime 7 case did not agree as favorably 
and at present no explanation is available for this problem. 

The drag increase was much higher for the glaze cases when com- 
pared to those for the rime shapes, however the effect of sur- 
face roughness did not appear to be as crucial. 

For the climb cases at the lower lift coefficients, the 
rime case shows a more negative pitching moment while at the 
higher C^'s the crend was reversed. For the glaze cases, however, 
more positive Cm's were found at the lower lift coefficients and 
more negative at higher C^'s. The change in C m is almost negli- 
gible for the cruise case at low C^'s. 

A run summary has been included at the end of this report. 

It is arranged by configuration; 1) clean, 2) rime 3, 3) glaze 3, 
*1) rime 7, and 5) glaze 7. These tables are further subdivided 
into rough and smooth. Following these are the cumulative 
vs. a, C(j vs. C£ and vs. C£ plots. Lastly, the C p distribu- 
tions are included ar.d are organized in the same order as the 
run summary tables. 
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An experimental program was conducted in which icc accre- 
tions on airfoils were measured, rnd simulated shapes tested, to 
determine tno degradation in airfoil section performance. The 
research was conducte \ in the NASA Icing Research Tunnel on a 
full-scale general aviation airfoil of construction and surface 
conditions tvpical of that found in industry. The airfoil model 
and simulated ice accretions were instrumented to obtain surface 
pressures which provided additional Insight into the flowfield 
behavior. 

As a result of this work some important observations can be 

made: 

1) The model must span the entire tunnel to obtain 
accurate 2-D results. Data accuracy could be 
improved by increased accuracy in the angle of 
attack measurement. 

2) The approximate ice shape modelling technique of 
a smooth shape with roughness added accurately 
simulated the ice in most casesl*2. The addition 
of surface roughness was particularly important 
in the rime ice simulation. 

3) Measured C p distributions on the ice shapes showed 
extremely severe adverse gradients. 

4) The data showed a significant increase in drag 
for all ice shapes tested. However, the maximum 
lift coefficient actually increased in one case 
and in all cases was higher with the simulated 
ice than expected. A zero lift angle of attack 
shift was also noted. 

This test demonstrated the capability of the IRT and the 
present ice simulation methods to generate aerodynamic data. 
While the comparisons of the aerodynamic performance of simu- 
lated and actual ice in the IRT appear good, an external check 
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on the quality of these data are still needed. The results pre- 
sented here represent a first step in generating the needed data 
base on the effect of ice accretion on airfoil performance. 
Combined with future results, these data will provide a means 
by which current analysis methods can be validated, as well as 
provide much needed guidance in the development of new methods. 
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CLEAN CONFIGURATION 


RUN 

AOA 1 

VEL 2 

TEMP 3 

PALT 4 

Cf 

Cd 

Cm 

2 

2.6 

99.25 

58.0 

1297.3 

0.6288 

0.0112 

-0.0848 

3 

:.6 

99.50 

60.0 

1304.4 

0.6540 

0.0126 

-0.0868 

4 

-0.4 

96.55 

60.0 

1278.0 

0.2677 

0.0116 

-0.0726 

5 

-3.4 

98.95 

49.0 

1295.2 

-0.0463 

0.0104 

-0.0675 

6 

1.6 

99.90 

58.0 

1303.4 

0.4313 

0.0120 

-0.0667 

7 

3.6 

98.75 

62.0 

1285.9 

0.7572 

0.0126 

-0.0909 

8 

5.6 

103.25 

51.0 

1319.2 

0.9228 

0.0139 

-0.0819 

9 

7.6 

99.80 

55.0 

1289.3 

1.1485 

0.0149 

-0.0855 

10 

8.6 

97.35 

49.0 

1260.6 

1.1701 

0.0156 

-0.0778 

11 

9.6 

100.00 

58.0 

1280.1 

1.1889 

0.0165 

-0.0711 

12 

10.6 

98.85 

52.0 

1269.1 

1.2751 

0.0180 

-0.0783 

13 

11.6 

98.75 

59.0 

1263,4 

1.3638 

0.0183 

-0.0942 

14 

2.6 

98.15 

79.0 

1179.1 

0.5894 

0.0130 

-0 . 0707 

15 

7.6 

102.20 

71.0 

1218.6 

1.0544 

0.0177 

-0.0507 

16 

-0.4 

100.15 

34.0 

1202.3 

0.2900 

0.0113 

-0.0685 

17 

11.6 

99.50 

49.0 

1196.9 

1.3699 

0.0250 

-0.0915 

18 

12.6 

98.95 

53.0 

1188.1 

1.4232 

0.0265 

-0.1044 

19 

13.6 

98.05 

51.0 

1170.2 

1.3910 

0.0344 

-0.0958 

20 

14.6 

100.70 

52.0 

1194.5 

1.3963 

0.0387 

-0.1239 

21 

2.6 

108.25 

53.0 

1275.6 

0.6606 

0.0118 

-0.0746 

22 

2.6 

122.15 

54.0 

1429.1 

0.6312 

0.0122 

-0.0673 

56 

2.6 

100.05 

60.0 

1232.5 

0.6609 

0.0122 

-0.0758 

57 

7.6 

99.15 

58.0 

1224.3 

1.2974 

0.0190 

-0.0809 

58 

-3.4 

81.90 

56.0 

1079.4 

0.0086 

0.0122 

-0.0835 

59 

-3.4 

82.45 

54.0 

1074.1 

0.0015 

0.0114 

-0.0822 

60 

-0.4 


54.0 


0.3447 

0.0113 

-0.0864 

61 

2.6 

84.30 

54.0 

1089.5 

0.6363 

0.0121 

-0.0751 

62 

4.6 

82.40 

54.0 

1073.6 

0.9467 

0.0129 

-0.0979 

63 

6.6 

80.35 

54.0 

1059.3 

1.0343 

0.0166 

-0.0668 

64 

8.6 

83.35 

57.0 

1076.8 

1.1823 

0 . 0218 

-0.0576 

65 

10.6 

82.15 

57.0 

1077.3 

1.3309 

0.0263 

-0.0595 

66 

11.6 

81.50 

57.0 

1078.7 

1.4416 

0.0324 

-0.0529 


^De^rees 

^Knots 

3op 

4 Ft 
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RIME 3 SMOOTH CONFIGURATION 


RUN 

AOA 

VEL 

TEMP 

PALT 

c* 

c d 

C™ 

114 

-3.4 

111.65 

45.0 

1362.1 

-0.1004 

0.0233 

-0.0851 

115 

-0.4 

112.90 

53.0 

1373.7 

0.2792 

0.0131 

-0.0862 

116 

2.6 

113.40 

50.0 

1387.5 

0.6182 

0.0122 

-0.0669 

117 

4.6 

113.25 

50.0 

1380.8 

0.8742 

0.0141 

-0.0725 

118 

6.6 

113.40 

48.0 

1380.3 

1.1129 

0.0179 

-0.0630 

119 

8.6 

111.35 

49.0 

1499.4 

1.2627 

0.0220 

-0.0439 

120 

9.6 

113.65 

48.0 

1512.4 

1.4184 

0.0252 

-0.0497 

121 

10.6 

112.00 

48.0 

1491.7 

1.5428 

0.0290 

-0.0534 

122 

11.6 

112.65 

49.0 

1495.1 

1.5689 

0.0347 

-0.0386 

123 

12.6 

112.65 

50.0 

1491.1 

1.6200 


-0.0453 

124 

13.6 

111.25 

50.0 

1477.0 

1.6669 


-0.0598 


RIME 3 ROUGH CONFIGURATION 


91 

-3.4 

113.75 

63.0 

1448.1 

-0.0705 

0.0210 

-0.0875 

92 

-0.4 

114.25 

51.0 

1451.1 

0.2780 

0.0149 

-0.0765 

93 

2.6 

113.40 

49.0 

1437.2 

0.6161 

0.0156 

-0.G656 

94 

-0.4 

113.10 

47.0 

1379.6 

0.3008 

0.0148 

-0.0800 

95 

4.6 

112.10 

48.0 

1355.2 

0.9274 

0.0190 

-0.0673 

96 

6.6 

111.60 

49.0 

1348.6 

1.0707 

0.0265 

-0.0578 

97 

8.6 

111.90 

48.0 

1337.6 

1.1951 

0.0348 

-0.0372 

98 

9.6 

112.65 

48.0 

1343.8 

1.2996 

0.0432 

-0.0419 

99 

10.6 

115.10 

48.0 

1376.6 

1.3705 

0.0517 

-0.0328 

100 

11.6 

114.00 

48.0 

1366.1 

1.4113 

0.0714 

-0.0380 

101 

12.6 

112.75 

48.0 

1346.3 

1.4463 


-0.0315 

102 

13.6 

112.30 

49.0 

1340.4 

1.4043 


-0.0458 
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RUN 

AOA 

VEL 

TEMP 

PALT 

c l 

Cd 

c m 

m 

37 

2.6 

107.60 

74.0 

1297.5 

0.6681 

0.0190 

-0.0738 

38 

-0.4 

106.90 

57.0 

1286.9 

0.3148 

0.0157 

-0.0738 

39 

7.6 

108.00 

57.0 

1299.6 

1.0679 


-0.0287 

40 

7.6 

108.50 

57.0 

1307.6 

1.1208 


-0.0383 

41 

7.6 

109.00 

57.0 

1307.4 

1.0983 

0.0668 

-0.0430 

42 

8.6 

108.15 

58.0 

1309.5 

1.2054 


-0.0301 

43 

10.6 

108.70 

57.0 

1314.7 

1.1935 


-0.0356 

44 

11.6 

109.65 

57.0 

1321.6 

1.1526 


-0.0433 

45 

5 . C 

108.15 

55.0 

1312.4 

1.0114 

0.0385 

-0.0487 

46 

-3.4 

107.65 

55.0 

1300.7 

-0.0427 

0.0248 

-0.0911 


GLAZE 3 ROUGH CONFIGURATION 


23 

2.6 

108.45 

62.0 

1283.5 

0.5998 

0.0161 

-0.0653 

24 

2.6 

108.80 

62.0 

1287.5 

0.5993 

0.0180 

-0.0650 

25 

-3.4 

108.45 

59.0 

1283.5 

-0.0492 


-0.0798 

26 

-3.4 

108.15 

56.0 

1280.0 

-0.0445 

0.0247 

-0.0907 

27 

-0.4 

107.90 

58.0 

1277.6 

0.2617 

0.0177 

-0.0762 

28 

4.6 

108.25 

58.0 

1278.1 

0.8034 

0.0271 

-0.0530 

29 

5.6 

108.85 

62.0 

1274.7 

0.9575 

0.0343 

-0.0536 

30 

6.6 

107.55 

57.0 

1254.0 

1.0038 


-0.0357 

31 

7.6 

107.20 

57.0 

1238.2 

1.1051 

— 

-0.0371 

32 

8.6 

108.70 

55.0 

1247.2 

1.1070 


-0.0213 
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FIGURE 1C. G3 ICE SHAPE 
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